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METHOD AND APPARATUS FOR DETERMINING AN EMG-SIGNAL 

5 The present invention is directed to a method according to the preamble to claim 1, 
a device according to the preamble to claim 7 and a computer program product ac- 
cording to the preamble to claim 13. 

Sensing of the EMG-activity in a patient's diaphragm by placing a catheter with a 
10 number of electrodes down the esophagus is a known technique, which is described 
in, among others, US 5,671,752. The EMG-signals, which are received by handling 
the signals, can be used in connection with mechanical ventilation of patients, which 
is described in, among others, US 5,820,560 and WO 98/48877. 

1 5 Sensing of the EMG-activity from the diaphragm can even be done outside the body 
by placing electrodes on the patient, described, for example, in US 4,248,240. 

The physiological activity in the diaphragm generates a relatively weak electrical 
signal. This is so in particular if it is compared with the considerably stronger 

20 physiological (and electrical) activity in the heart (EKG-activity). Therefore, the 
present desire is to attain in the best way the highest quality possible for the signal 
handling of the raw signal which the sensors detect, so that the resulting EMG- 
signal in the highest degree possible corresponds to the physiological activity. 
This is evident even in WO 01/03579. In WO 01/03579 it is assumed that the elec- 

25 trodes location in relation to the center of the diaphragm is known. Then the elec- 
trodes are measured based on location and symmetry, in which the EKG signal is 
measured in a traditional way. 

Previously known methods for compensating for the EKG signal include, among 
30 others, using a band-pass filter which filters out the frequencies where the EKG sig- 
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nal normally appears. It is also known to measure the EKG signal separately and 
then remove an equivalent signal from the measured EMG signal. 

A specific problem which is present with interference from EGK-activity is that the 
5 frequency spectrum for the physiological EKG-activity partially overlaps the fre- 
quency spectrum for the physiological EMG-activity, Moreover, there is the fact that 
the EKG-disturbances appear at different time points in relation to the breathing cy- 
cle. 

10 None of the above-described methods consider the actual disturbance the EKG sig- 
nal creates in a particular measuring situation in a particular patient. This distur- 
bance also varies with time. 

The known techniques also fail to consider that the respective frequency spectrum 
15 for EMG and EKG can vary between measurements in different people. 

It is a goal with the present invention to provide a device and a method which lead 
to an improvement of known methods and procedures to filter out EMG-signals. 

20 Another goal with the invention is to provide a device which can lead to an im- 
proved filtering of a raw signal to receive an EMG-signal. 

This goal is reached in combination with the invention in that the method above 
comprises the characteristics in the characterizing part of claim 1, and the device 
25 further includes the characteristics which are specified in the characterizing part of 
claim 6. 

The basis of the invention is that an EKG-signal and an EMG-signal are first esti- 
mated out of the raw signal. An appropriate detection window for determining the 
30 EMG-signal can then be determined from the estimated signals. In this connection. 
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the detection window can be detennined with reference to the sensitivity, width 
(frequency) or position in the frequency plane. The determination can even be made 
from a combination of two or more of these different possibilities. 

5 In that way the invention makes possible an adaptable filtering of the EMG signal 
by measuring the actual frequency range for the respective EMG and EKG signal at 
each opportunity when the border frequency for the filter is to be detennined. 

An example of determining sensitivity is to first identify the frequency, where the 
1 0 EKG-interference becomes weaker than the EMG-signal. For that overlapping fre- 
quency region the EMG-signal is determined according to a first criterion (harder 
filtering, subtraction of estimated EKG-interference, etc.). In the non-overiapping 
part, the EMG-signal is determined according to another criterion. In other words, 
the EMG- window is divided in two (or more) sub-windows, each with different sig- 
1 5 na! handling conditions. 

An example of variation of the width of the detection window is to first identify the 
frequency where the EKG-disturbance becomes weaker than the EMG-signal. The 
size of the detection window is then chosen to correspond to the distance from this 
20 frequency and the remaining part of the detected EMG-spectrum. 

An example of variation of the position of the detection window is to first identify 
the frequency where the EKG-disturbance becomes weaker than the EMG-signal. 
Then the entire detection window is displaced so that this starting point coincides 
25 with the determined frequency. 

In a further development of the method, consideration is taken to the noise level. In 
this connection, the detection window can also be limited to fall within the region 
where the EMG-signal is stronger than the noise. 



30 
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An example of the method and device according to the invention is described in 
connection with the figures. 



Fig. 1 shows a device according to the invention connected to a patient who is re- 
5 ceiving respiratory therapy. 

Fig. 2 shows a first example of an electrode connection to receive a raw signal. 
Fig. 3 shows a second example of an electrode connection to receive a raw signal. 
Fig. 4 shows schematically how the EKG-activity can affect detection of the EMG- 
activity over time, 

10 Fig. 5 shows schematically how frequency spectra for EKG and EMG partially over- 
lap and how the method according to the invention can be designed, and 
Fig. 6 shows a flow chart for an analysis unit according to the invention. 

Fig. 1 shows a device 2 for determining EMG-signals according to the invention. 

15 The device 2 can in a known way be connected to a patient 4 via a catheter 6 with a 
number of electrodes 8A, 8B, 8C, 8D in the tip (four electrodes are shown, but the 
number can be larger or smaller). By placing the catheter 6 in esophagus (gullet) of 
the patient 4, the electrodes 8A, 8B, 8C, 8D can be placed in different locations in 
the diaphragm 10 (whose size is exaggerated in the figure to indicate the relative 

20 placement of the electrodes 8A, 8B, 8C, 8D). In an analysis unit 12 in the device 2, 
filtering and analysis of the raw signal from the catheter 6 is done to extract the 
highest quality EMG-signal possible, i.e. an EMG-signal which corresponds to the 
physiological EMG-activity. 

25 In this connection, the raw signal can be received in many different ways. Fig. 2 and 
Fig. 3 exemplify two ways. In Fig. 2 it is evident that the electrodes 8A, 8B, 8C, 8D 
can be connected in pairs via three couplers 14 A, 14B, 14C and thereby give rise to 
a three channel raw signal (with e.g. nine electrodes, eight channels are received in a 
corresponding way). 
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In Fig. 3 an example is shown where the respective electrodes 8A, 8B, 8C, 8D are 
connected to a reference electrode 16 (which e.g. can be grounded) via four couplers 
14D, 14E, 14F, 14G. This gives rise to a four channel raw signal (for eight channels 
in this configuration, eight electrodes and a reference are consequently required), 

5 

More information regarding the catheter, the sensors and the entire process to cap- 
ture raw signals from the diaphragm via the esophagus is given in e.g. US 5,671,752 
and WO 01/03579. As noted, electrodes connected outside the body can be used in- 
stead for receiving EMG-signals completely non-invasively. 

10 

The patient 4 can even be connected in a conventional way to a ventilator system 1 8, 
which can be connected to a device 2. The respiratory therapy, which is given via 
the ventilator system 1 8, can in that way be influenced by the EMG-signal which is 
extracted from the raw signal from the diaphragm 10. This influence can be done in 
15 many different ways, of which some are described in US 5,820,560 and 
WO 99/43374. 

The present invention is directed to a device 2 and, to be precise, the analysis unit 
12. The analysis unit 12 filters EMG-signals out of the raw signal from the catheter 
20 6. In this connection, a number of signal channels are used, as noted above. 

To receive the highest quality EMG-signals possible, the filtering in the analysis 
unit 12 is done according to the method described below, which can be performed 
analog, digital or in a combination of these and realized in hardware, software of a 
25 combination of these. 

One of the biggest interferences in measuring physiologic EMG-activity (determin- 
ing of a representative EMG-signal) comes from the heart. Fig. 4 shows schematic- 
ally how a number of EMG-activities 18 (representing the electrical activity of the 
30 diaphragm during inhalation) periodically are flooded by considerably stronger 
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EKG-activities 20. Since breathing (the EMG-activity 1 8) as well as heart beat (the 
EKG-activity 20) have different periodicity (which in itself is variable), the signals 
do not coincide in a regular or predictable way, which makes the determination of 
the EMG-signal out of the raw signal more difficult. 

5 

Fig. 4 also shows a noise 22, which also contributes to making the measurement of 
EMG-activity 18 more difficult. 

Fig. 5 shows schematically how an EMG-signal 24 and an EKG-signal 26 partially 
10 contain overlapping frequency regions in their respective (physiological) frequency 
spectra. With a frequency analysis limited to the frequency region of the EMG- 
signal 24, it cannot therefore be excluded that the EKG-signal 26 gives a certain 
contribution. At the same time, a limiting of the frequency region which is analyzed 
should result in a larger uncertainty in the filtering of the EMG-signal 24 as the total 
15 EMG-energy risks being too small to effectively be detected. 

Fig. 5 shows also a noise 28, which in this case comprises a white background noise 
with relatively constant energy content within the entire frequency region. Like the 
EKG-signal 26, the noise 28 results in the frequency region for detection of the 
20 EMG-signal 24 being limited. 

With help from Fig. 4 and Fig. 5, different examples of the method according to the 
invention can be described. As is evident from Fig. 4, during each inhalation 18, one 
or a number of EKG-activities 20 regularly overlap. By analyzing the EKG-signal 
25 26 within a frequency interval feKo and at the same time estimating the frequency 
spectrum of the EMG-signal 24, the EKG-signaPs 26 appearance within the over- 
lapping part of the EMG-signal's 24 frequency spectrum can also be estimated. This 
in turn can be used to influence how the EMG-analysis is to advance. 
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A first way is then to determine a lower frequency border fcl, where the EKG- 
signal 26 and the EMG-signal 24 intersect, and let this frequency fcl constitute the 
lower border in an analysis window for the EMG-signal (EMG-window). The upper 
border for the EMG-window can be determined in a corresponding way to an upper 
5 border frequency fc2, where the EMG-signal 24 intersects the noise signal 28. The 
lower border frequency fcl is set in the same way in the absence of EKG signal. 

According to this process, the EMG-window should be adjusted to its frequency 
width for, in principle, each detection. Yet this analysis involves at the same time 
10 that the noise energy contribution varies since the EMG-window width varies. As- 
suming that the noise is relatively constant withm the frequency spectrum where the 
EMG-window width is going to vary, that varying contribution from the noise can 
be compensated for each determination. 

15 It can be noted that it is not necessary to directly use the lower border frequency fcl 
as a starting point for the EMG-window. The determination of the EKG-signal 26 
also gives a determination regarding the slope for the EKG-signal 26 at the lower 
border frequency fcl. In this connection, if the EKG-signal 26 is flat, there can be 
reason to displace the EMG-window even more as the contribution from a flat 

20 EKG-curve decreases more slowly than the contribution from a steep EKG-curve. In 
other words, the EMG window is displaced nearer the EMG signal's middle fre- 
quency fc so that the contribution from the EKG signal reduces to increase the 
safety margin. 

25 An alternative way to analyze the EMG-signal according to the invention is to in- 
stead use an EMG-window with constant width, but use a lower border frequency 
fcl as a starting point for the EMG-window (which then can extend longer than to 
the upper border frequency fc2). In that way, a substantially constant contribution 
from the noise-signal 28 will be received. 

30 
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In the same way as above, the slope (derivative) of the EKG-curve 26 at the lower 
border frequency fcl can also be used to further finely adjust the placement of the 
. EMG-window. 

5 Yet another alternative way to analyze the EMG-signal 24 according to the inven- 
tion is to use a constant EMG-window, which is big enough to contain the entire 
EMG-signal, but which falls outside the interval fEKO- From the estimate of the 
EKG-signal 26 which is made, the contribution which the EKG-signal entails for 
each frequency (or total) can then be determined, which then is subtracted from the 
10 EMG-window. In principle, an "undisturbed" EMG-signal can be determined in this 
way (after subtraction of the contribution from the noise). 

It is also possible to use the estimate of the EKG-signal 26 in such a way that the 
EMG-window is divided into sub-windows with different sensitivity or weighing. 

15 

That which is mentioned above can also be performed for the entire noise signal 28 
(which can be estimated within a frequency region over the highest frequencies of 
the EMG-signal). 

20 There are further ways to use the determination. In principle, the above has an im- 
mediate adaptation of the EMG-window occurred (in real time or with backlog) on 
the raw signal, i.e. a use of each specific EMG-activity (breath). 

A slower adjustment of the EMG-window is possible, too, where trends over time 
25 are used to successively adjust the EMG-window. An advantage of slower adapta- 
tion is that patient-specific properties in the raw signal can be used in another way 
than with the immediate adjustment. 

It can even be noted that the frequency spectra for the EMG-signal 24 as well as the 
30 EKG-signal 26 are not constant but can vary (within physiological borders). 
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In particular, the middle frequency for the frequency spectrum of the EMG-signal 
24 can be displaced. This principally occurs dependent on changes within the pa- 
tient, e.g. muscle fatigue or a reaction to the amount of relaxing medications which 
5 may have been given. 

Minding the middle frequency can therefore be used in many different ways for the 
treatment. For example, an alarm can be generated at too high or fast variations in 
displacement of the middle frequency. Increased breathing support can be suggested 
10 to the operator or initiated automatically if the middle frequency is displaced beyond 
a border limit. The level of relaxing medication can be watched by studying the 
middle frequency variation (and even be used to maintain a certain level of calm, i.e. 
for regulating the dosage of the calming medicine). 

15 Figure 6 shows a logical flow chart of an analysis unit 12 adapted to perform the 

method according to the above. The raw signal from the electrodes comes in via the 
catheter 6 to the analysis unit 12. After an AC-connection unit 30 the signal passes 
to an adaptive band-pass filter 32. The signal out from the band-pass filter 32 is a 
raw signal where the EKG signal is filtered out. This signal goes fiirther to energy 

20 calculation which gives the desired EMG signal. The band-pass filter 32 uses a filter 
coefficient 0, which is measured according to the description below. The signal 
from the AC-connection unit 30 and the output signal from the band-pass filter 32 
are even fed to a time signalling block 34 whose output signal is a logical signal 
with four possible values: "not valid", "EMG", "EKG" and "noise", "not valid" 

25 means that no estimate of the signal will be made (if, for example, EKG is synchro- 
nous with EMG). If the value is EMG, EKG or "noise" it means that the respective 
spectrum can be estimated. 

For estimating the different signals the output signal is fed from the AC-connection 
30 unit 30 to three different measuring units. A first 36 of these measuring units calcu- 
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lates the sum of noise and the EKG signal. A second measuring unit 38 calculates 
the sum of noise and the EMG signal. A third measuring unit 40 calculates noise 
only. The three measuring units are controlled by the logical ouQ>ut signal from the 
time signalling block 34. An output signal from the third measuring unit 40 is used 
5 for calculating energy. The output signal from the third is used to remove the en- 
ergy level equivalent to the noise signal. 

The ou^ut signal from the second measuring unit 38 and the third measuring unit 
40 can even be used to calculate the middle frequency for EMG, which in that case 

10 is done in a calculating unit 44. For example, the middle frequency can be deter- 
mined from the point of balance /mean for the filtered spectrum or alternatively 
maximum amplitude. As complimentary information the range of measurements 
/variance of the EMG signals spectra can be used. The variations of the middle fre- 
quency is a measure of muscle fatigue, and can also indicate that the patient is se- 

15 dated, and in that case, how much. If the patient becomes exhausted an alarm can be 
sent, for example, or possibly ventilation assistance can increase. 

The output signal from all the measuring blocks 36, 38 and 40 is fed to a filter opti- 
mising block 42 which calculates the filter coefficient 0 for the band-pass filter 32. 

20 

Alternative methods exist to measure the spectrum in the blocks 36, 38 and 40. One 
alternative is to use Fourier analysis. Another alternative is to make an identifica- 
tion/adjustment of a model of the signal, for example, with the recursive method of 
least squares with forgetting factor, where the model's parameters are continually 
25 adapted to changes in the time signal. 

Depending primarily on the cutting between the EMG spectra and the EKG spectra 
the filter order can be chosen. The filter should not be chosen unnecessarily sharp 
since this impacts the energy calculation in the time signal. Preferably, a steep filter 
30 with high ranking is chosen if the energy in the EKG signal decreases gradually in 
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the region around the cut at the same time that the energy in the EMG signal in- 
creases gradually (that is, if the curves are flat). If the respective energy decreases or 
increases faster, a filter with lower ranking can be chosen. 

In a later step the signal energy is determined. Since the bandwidth varies even the 
noise's energy contribution will vary. For a chosen pass band the energy is calcu- 
lated for noise plus EMG. Therefore the noise's energy spectrum is used to compen- 
sate so that a constant energy contribution is received which corresponds to a con- 
stant pass band. 

Alternatively one could use a constant width in the pass band, since the critical re- 
gion is that around the lower breaking frequency. 

The skilled worker naturally understands that the flow chart shovm in figure 6 is in- 
tended as a greatly simplified illustration of a working example and that an analysis 
unit to perform the method can be achieved in many different ways. 

Preferably, the analysis unit comprises a computer program which directs the input- 
ting units' function. 

A number of algorithms and filter categories to use for estimations and determina- 
tions according to the above are known to a person skilled in the art and can be 
combined in many obvious ways. There is therefore no reason to describe them in 
more detail in this context. 



